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WHAT IS CLAIMED IS : 

1 1 . A method for maximum a posteriori (MAP) decoding of an input 

2 information sequence based on a first information sequence received through a channel, 

3 comprising: 

4 iteratively generating a sequence of one or more decode results starting with an 

5 initial decode result; and 

6 outputting one of adjacent decode results as a decode of the input information 

7 sequence if the adjacent decode results are within a compare threshold. 

1 2. The method of claim 1, wherein the iteratively generating comprises: 

2 a. generating the initial decode result as a first decode result; 

3 b. generating a second decode result based on the first decode result and a model 

4 of the channel; 

5 c. comparing the first and second decode results; 

6 d. replacing the first decode result with the second decode result; and 

7 e. repeating b-d if the first and second decode results are not within the compare 

8 threshold. 

1 3. The method of claim 2, wherein the generating a second decode result 

2 comprises searching for a second information sequence that maximizes a value of an 

3 auxiliary function. 

1 4. The method of claim 3, wherein the auxiliary function is based on the 

2 expectation maximization (EM) algorithm. 

1 5. The method of claim 4 ? wherein the model of the channel is a Hidden 

2 Markov Model (HMM) having an initial state probability vector 71 and probability density 

3 matrix (PDM) of P(X,Y), where XeX, YeY and elements of P(X 5 Y) 5 

4 py(X 5 Y) = Pr(j,X 5 Y | i), are conditional probability density functions of an information 

5 element X of the second information sequence that corresponds to a received element Y 

6 of the first information sequence after the HMM transfers from a state i to a state j, the 

7 auxiliary function being expressed as: 
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8 Q(X^,X^ p ) = ^^(z,X^ p) Y, T )log(T(z,X^ Y, T )) 5 where p is a number of 

T 

9 iterations,^, Xj , Y?) = 71 ^f^Pw ,( X t> Y t) > T is a number of information elements 

10 in a particular information sequence, z is a HMM state sequence ij , ic s is the probability 

11 of an initial state io 5 Xj is the second information sequence, is a second information 

12 sequence estimate corresponding to a pth iteration, and Y, T is the first information 

13 sequence. 

1 6. The method of claim 5, wherein the auxiliary function is expanded to be: 

2 Q( X i. X i) = EZI YUj(X;, p )log(p ij (X t ,Y t )) + C 

t=l i=l j=l 

3 where C does not depend on X[ and 

4 T^( X np)= a .( x r [ » l ' Y . , " l )Pij(X,,p,Y t )(3j(X^ p ,Y t T +l ) 

5 where a i (Xj p , Y, T ) and Pj(X^ +I p , Y,+, ) are the elements of forward and backward 

6 probability vectors defined as 

7 a(X;, Y, , ) = 7 l flP(X i ,Y i ),andp(X:,Y 1 T ) = nP(X j ,Y j )l, n is an 

i=l j=t 

8 initial probability vector, 1 is the column vector of ones. 

1 7. The method of claim 6, wherein a source of an encoded sequence is a 

2 trellis code modulator (TCM), the TCM receiving a source information sequence l[ and 

3 outputting Xf as an encoded information sequence that is transmitted, the TCM defining 

4 X t = g t (S t ,I t ) where X t and I t are the elements of Xf and if for each time t, respectively, S t 

5 is a state of the TCM at t, and g t (.) is a function relating X t , to I t and S t , the method 

6 comprising: 

7 generating, for iteration p+1, a source information sequence estimate ij p+1 that 

8 corresponds to a sequence of TCM state transitions that has a longest cumulative distance 

9 L(Sm) at t = 1 or L(So), wherein a distance for each of the TCM state transitions is 
10 defined by L(S t ) = L(S t +i) + m( I t+1 (S t +i))for the TCM state transitions at each t for 
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11 t = 1, T and the cumulative distance is the sum of m(I ( (S t )) for all t 5 m(I t (S t )) being 

12 defined as 

13 m(l t (S t ))=X X Yg j fc)logp^(YjX t (S t )),foreacht=l,2,...,T > where 

14 X t (S t ) = gt(S t? I { (S t )), n c is a number of states in an HMM of the channel and 

1 5 Pc,ij(Y t I X t (S t )) are channel conditional probability density functions of Y t when X t (S t ) is 

16 transmitted by the TCM, if p+1 being set to a sequence of I t for all t. 

1 8. The method of claim 7, wherein for each t = 1,2, T, the method 

2 comprises: 

3 generating m( I t (S t )) for each possible state transition of the TCM; 

4 selecting state trajectories that correspond to largest 

5 L(S { ) = L(S t +i) + m( I t+1 (S t +0) for each state as survivor state trajectories; and 

6 selecting I t (S t )s that correspond to the selected state trajectories as I t , p +i(S t ). 

1 9. The method of claim 8, further comprising: 

2 a. assigning L(S T )=0 for all states at t = T; 

3 b. generating m( I t (S t )) for all state transitions between states S t and all possible 

4 states S t+ i; 

5 c. selecting state transitions between the states S t and S t +i that have a largest 

6 L(S t ) = L(S t+ i ) + m(I t+1 (S t +0) ai *d I t+ , (S t +0 that correspond to the selected state 

7 transitions; 

8 d. updating the survivor state trajectories at states S t by adding the selected state 

9 transitions to the corresponding survivor state trajectories at state S t +i ; 

10 e. decrementing t by 1 ; 

1 1 f. repeating b-e until t = 0; and 

12 g. selecting all the I t (S t ) that correspond to a survivor state trajectory that 

1 3 corresponding to a largest L(S t ) at t = 0 as ij x . 

1 10. The method of claim 6, wherein the channel is modeled as P C (Y | X) = 

2 P C B C (Y | X) where P c is a channel state transition probability matrix and B C (Y | X) is a 



A 
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3 diagonal matrix of state output probabilities, the method comprising for each t = 1, 2 5 

4 T: 

5 generating y tJ ( I J, ) = ai ( Y/ I Ij, )P,( Y t T +1 | p ); 

6 selecting an I t (S t ) that maximizes L(S t )=L(S t +i )+ m( I t+1 (S t +i)), where 

7 m( 1 1 (S t )) is defined as 

8 m( 1 1 (S t )) = Yj ft i( l l? )Pj(Y t | X t (St)) ? n c being a number of states in an HMM of 

i=l 

9 the channel; 

10 selecting state transitions between states S t and S t +i that corresponds to a largest 

11 L(S t ) = L(S t+1 ) + m(i t+1 (S t+1 ));and 

12 forming survivor state trajectories by connecting selected state transitions. 

1 11. The method of claim 1 0, further comprising: 

2 selecting I t (S t ) that corresponds to a survivor state trajectory at t = 0 that has the 

3 largest L(S t ) as I 1 r p+1 for each pth iteration; 

4 comparing lj p andl^ p+l ; and 

5 outputting I* p+] as the second decode result if I^ p andl* , are within the 

6 compare threshold. 

1 1 2. A maximum a posteriori (MAP) decoder that decodes a transmitted 

2 information sequence using a received information sequence received through a channel, 

3 comprising: 

4 a memory; and 

5 a controller coupled to the memory, the controller iteratively generating a 

6 sequence of one or more decode results starting with an initial decode result, and 

7 outputting one of adjacent decode results as a decode of the input information sequence if 

8 the adjacent decode results are within a compare threshold. 

1 13. The decoder of claim 12, wherein the controller: 

2 a. generates the initial decode result as a first decode result; 

3 b. generates a second decode result based on the first decode result and a model 

4 of the channel; 
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5 c. compares the first and second decode results; 

6 d. replaces the first decode result with the second decode result; and 

7 e. repeats b-d until the first and second decode result are not within the compare 

8 threshold. 

1 14. The decoder of claim 13, wherein the controller searches for information 

2 sequence that maximizes a value of an auxiliary function. 

1 15. The decoder of claim 14, wherein the auxiliary function is based on 

2 expectation maximization (EM). 

1 16. The decoder of claim 15, wherein the model of the channel is a Hidden 

2 Markov Model (HMM) having an initial state probability vector n and probability density 

3 matrix (PDM) of P(X,Y), where XeX, YeY and elements of P(X,Y), 

4 Py(X,Y) = Pr(j 3 X,Y | i), are conditional probability density functions of an information 

5 element X of the second information sequence that corresponds to a received element Y 

6 of the first information sequence after the HMM transfers from a state i to a state j, the 

7 auxiliary function being expressed as: 

8 Q(Xj , X J p ) = £ ¥(z, X^ , Y, T ) log(y (z, Xf , Y, T )), where p is a number of 

T 

9 iterations, ^(z, Xf, Y, T ) = n io fjp it li( (X t ,Y t ) , T is a number of information elements 

10 in a particular information sequence, z is a HMM state sequence ij , % . is the probability 

11 of an initial state io, X^ is the second information sequence, X^ p is a second information 

12 sequence estimate corresponding to a pth iteration, and Y, T is the first information 

13 sequence. 

1 17. The decoder of claim 16, wherein the auxiliary function is expanded to be: 

2 Q(X^,X^ p ) = i t Z Y«jj(X^) log ( Pij (X t , YO) + C 

t=l i=l j=l 

3 where C does not depend on X, and 

4 Y (X^ p )=a, (XI; 1 , Y/" 1 ) Pij (X tjp , Y t ) Pj( X^ , Y t T +1 ) 

5 where a ■ (X J p , Y, T ) and Pj( X[ +1 p , Y t +, ) are the elements of forward and backward 

6 probability vectors defined as 
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7 a( XJ , Y; ) = 7t fl P(X f , Y s ) , and (3 (X* , , Y, T +1 ) = fl P(X } , Y } ) 1, n is an initial 

i=l j=t+l 

8 probability vector, 1 is the column vector of ones. 

1 18. The decoder of claim 17, wherein a source of an encoded sequence is a 

2 trellis code modulator (TCM), the TCM receiving a source information sequence ij and 

3 outputting Xj as an encoded information sequence that is transmitted, the TCM defining 

4 X t = g t (S t JO where X t and I t are the elements of and l[ for each time t, respectively, S t 

5 is a state of the TCM at t, and g t (.) is a function relating X t , to I t and S t , the controller 

6 generates, for iteration p+1, an input information sequence estimate ij p+] that 

7 corresponds to a sequence of TCM state transitions that has a longest cumulative distance 

8 L(Sm) at t = 1 or L(So), wherein a distance for each of the TCM state transitions is 

9 defined by L(S t +i) = L(S t +i) + m(I l+! (S t+ i)) for the TCM state transitions at each t for 

10 t = 1, T and the cumulative distance is the sum of m(I t (S t )) for all t, m(I t (S t )) being 

1 1 defined as 

12 m(l t (S t ))=f; £ Y tjij (l^ p )logp cij (Y t |X t (S t )),foreacht=l, 2, T, where 

13 X t (S t ) = g t (S t , I t (S t )) 5 n c is a number of states in an HMM of the channel and 

14 Pc,ij(Y t I X t (S t )) are channel conditional probability density functions of Y t when X t (S t ) is 

15 transmitted by the TCM, I^ p+I being set to a sequence of I ( for all t. 

1 19. The decoder of claim 1 8, wherein for each t = 1 , 2, T, the controller 

2 generating m( I t (S t )) for each possible state transition of the TCM, selecting state 

3 trajectories that correspond to largest L(S t ) = L(S t +i) + m(I t+I (S t +i)) for each state as 

4 survivor state trajectories, and selecting I t+I (S t +i)s that correspond to the selected state 

5 trajectories as I t +i 5 p+i(S t +i). 

1 20. The decoder of claim 19, wherein the controller: 

2 a. assigns L(St)=0 for all states at t = T; 

3 b. generates m(I t (S t )) for all state transitions between states S t and all possible 

4 states S t +i; 
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5 c. selects state transitions between the states S t and S t +i that have a largest 

6 L(S t ) = L(S t +i ) + m( I t+1 (S t +i)) and I t+1 (S t +i) that correspond to the selected state 

7 transitions; 

8 d. updates the survivor state trajectories at states S t by adding the selected state 

9 transitions to the corresponding survivor state trajectories at state S t +i; 

10 e. decrements t by 1 ; 

1 1 f. repeats b-e until t = 0; and 

12 g. selects all the I t (S t ) that correspond to a survivor state trajectory that 

13 corresponding to a largest L(S t ) at t = 0 as ij t . 

1 21. The decoder of claim 20, wherein the channel is modeled as P C (Y | X) = 

2 P C B C (Y | X) where P c is a channel state transition probability matrix and B C (Y | X) is a 

3 diagonal matrix of state output probabilities, for each t = 1, 2, T, the controller: 

4 generates YtiO = otj( Y/ | I' p )Pi( Y t T +1 | £ !>p ); 

5 selects an I t (S t ) that maximizes L(S t )=L(S t +i )+ m(I t+1 (S t +i)), where m(I t (S t )) 

6 is defined as 

7 m(i t (S t )) = J] Yt, i(ljp )Pj(Y t I X t (St)), ii c being a number of states in an HMM of 

i=I 

8 the channel; 

9 selects state transitions between states S t and S t +i that corresponds to a largest 

10 L(S t ) = L(S t+ , ) + m(I t+1 (Stf 0); and 

1 1 forms survivor state trajectories by connecting selected state transitions. 

1 22. The decoder of claim 2 1 , wherein the controller selects 1 1 (S t ) that 

2 corresponds to a survivor state trajectory at t = 0 that has the largest L(S t ) as ij } for each 

3 pth iteration, compares \J p andl^ p+1 , and outputs lj r p+l as the second decode result if 1^ 

4 andl^, are within the compare threshold. 



